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非 均 一 方差 的 GNSS 定 位 误差 模型 检验 - 
Reis ?, De, tock, PR, BA? 
(1. 中 国 科 学 院 国家 天 文 台 ， 北 京 100012; 2. 中 国 科学 院 大 学 ， 北 京 100049) 
摘要 : 在 全 球 导 航 卫 星系 统 (Global Navigation Satellite System, GNSS) 中 ， 定 位 
精度 是 一 个 重要 的 性 能 指标 ; 准确 的 定位 误差 模型 有 利于 理解 系统 的 工作 模式 , 同时 也 为 进 
一 步 提高 系统 性 能 提供 指导 。 当 卫星 的 测 距 误差 满足 相互 独立 和 方差 均一 的 条 件 时 , 定位 误 
差 等 于 测 距 误差 标准 差 乘 以 几何 精度 因子 .然而 实际 应 用 中 ,由 于 不 同 卫星 信号 的 传播 路 径 
的 不 同 , 测 距 误差 的 统计 特征 难以 相同 。 尤其 当 有 卫星 信号 受到 多 径 干扰 与 电离 层 闪 烁 影响 
时 ， 不 同 卫 星 间 的 测 距 误 差 会 有 显著 的 差异 。 本 文 介绍 了 一 种 利用 奇异 值 分 解 方法 建立 的 非 
均一 测 距 误 差 的 定位 误差 理论 模型 ; 利用 北京 站 与 香港 站 的 实测 数据 , 对 受到 多 径 效 应 与 电 
离 层 闪烁 影响 的 定位 误差 进行 分 析 。 结 果 表 明 其 统计 特性 与 非 均 一 方差 的 定位 误差 模型 给 出 
的 结果 一 致 ， 验 证 了 模型 的 正确 性 和 有 效 性 。 
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GNSS 卫星 定位 误差 分 析 中 ， 主 要 使 用 测 距 误差 均一 前 提 假 设 下 的 定位 误差 模型 ， 即 当 
卫星 的 测 距 误差 满足 相互 独立 和 方差 均一 的 同 分 布 条 件 时 , 定位 误差 等 于 测 距 误差 标准 差 乘 
以 几何 精度 因子 (Geometric Dilution of Precision, GDOP) 后。 几何 精度 因子 仅 与 卫星 
几何 分 布 情况 有 关 ， 它 与 用 户 接 收 机 与 可 视 卫星 形成 的 四 面体 或 多 面体 的 体积 成 反比 “。 有 
很 多 研究 工作 致力 于 参与 定位 选 星 工 作 ， 以 得 到 尽 可 能 小 的 GDOP 值 一 。 测 距 误差 是 对 卫星 
伪 距 测量 值 进行 误差 校正 后 得 到 的 伪 距 测量 误差 量 ,其 中 包含 未 剔除 干净 的 对 流 层 延迟 、 电 
离 层 延迟 、 多 径 干 扰 与 接收 机 噪声 等 误差 项 。 实 际 定 位 解 算 中 ， 由 于 卫星 间 的 相互 差异 、 传 
播 路 径 不 同等 因素 ， 不 同 卫星 测 距 误差 的 统计 特征 难以 相同 ， 无 法 满足 方差 均一 的 条 件 ”。 
尤其 当 某 颗 卫星 信号 受到 多 径 效应 或 电离 层 闪 烁 影响 时 ， 其 测 距 误 差 显著 增 大 “”。 围 绕 多 
径 干扰 、 电 离 层 内 烁 与 定位 误差 间 的 关系 展开 了 很 多 研究 工作 , 有 研究 给 出 统计 分 析 和 试验 
结果 ， 说 明了 多 径 干 扰 、 电 离 层 闪烁 对 定位 误差 存在 影响 “”。 但 尚未 给 出 统一 结论 描述 多 
径 效 应 、 电 离 层 闪 烁 与 定位 误差 之 间 的 关系 。 有 研究 提议 剔除 测 距 误差 大 的 星 ,但 这 样 做 会 
改变 卫星 星座 分 布 ， 可 能 会 恶化 几何 精度 因子 ”。 该 措施 对 定位 精度 的 维持 起 一 定 作用 ,但 
对 于 一 般 的 用 户 接收 机 并 不 实用 。 因 此 , 在 测 距 误 差 非 均 一 的 实际 定位 解 算 中 ,尤其 是 当 卫 
星 信号 受到 多 径 干 扰 或 电离 层 闪 烁 影响 时 ,如 何 保证 卫星 导航 定位 精度 , 如 何在 一 般 用 户 接 
收 机 中 实施 ， 需 要 考察 非 均一 测 距 误差 如 何 影响 定位 精度 。 已 有 研究 针对 GPS 定位 考察 了 非 
均一 伪 距 误差 对 定位 误差 的 影响 ， 给 出 了 定位 误差 模型 ”。 

本 文 旨 在 对 卫星 测 距 误差 非 均 一 时 的 定位 误差 模型 进行 验证 。 首 先 利 用 奇异 值 分 解 
(Singular Value Decomposition, SVD ) 方 法 对 几何 矩阵 进行 分 解 ， 推 导出 定位 误差 协 方差 
KAR, 给 出 定位 误差 模型 。 接 下 来 利用 北京 站 与 香港 站 的 实测 数据 ， 对 受到 多 径 效 应 与 电 
离 层 闪烁 影响 的 卫星 信号 进行 定位 解 算 , 将 定位 结果 与 非 均 一 测 距 误 差 下 的 理论 模型 进行 对 
比 和 分 析 。 最 后 给 出 了 结论 。 
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2 基于 SVD 方法 的 定位 误差 模型 


在 线性 代数 


有 数值 稳定 性 ， 相 比 易 受 舍 入 误差 影响 


P, SVD 是 分 解 和 矩阵 的 基本 方法 ， 利 用 
完全 等 效 的 。SVD 的 优势 在 于 它 能 够 给 出 矩阵 的 伪 逆 ， 从 而 推 
xk, M SVD 4 


SVD 求解 方程 在 理论 上 与 最 小 二 乘 是 
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8 解 的 简洁 、 自 约束 的 表达 


得 到 的 奇异 值 和 正 交 向 量 可 以 洞察 问题 的 几何 与 理论 内 涵 。 并 且 


问题 ， 它 是 最 小 二 乘 意义 下 的 最 优 近似 解 。 
GNSS 卫星 定位 中 ， 定 位 误差 方程 可 表示 为 : 


收 机 位 置 与 钟 差 的 未 知 修正 量 ，G 为 几何 矩阵 。 利 用 SVD 方 法 对 儿 何 矩阵 6 进行 分 解 ， 可 推 * 
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E 测 量 值 与 估 测 值 之 差 和 测 距 误差 ，[6x 6y 6z Sb]? 为} 
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定位 误差 协 方差 表达 式 ， 给 出 定位 误差 模型 。 


用 SVD 方 法 对 几何 矩阵 G 进 行 分 解 : 


G = Udiag(wi)VT 


6 矩阵 的 伪 逆 为 : 


其 中 
mii 


UAIVE I 


Gt = Vdiag(A,;)U™ 》 Ài = 1/%i 


Ui 1 U12 
Uz 1 Uz 2 


UNa UNz 


Vz 
Vy2 
V22 
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E 交 矩阵，U 的 列 向 量 和 V 的 列 向 量 分 另 


，SVD 求解 具 


方程 自身 奇异 的 最 小 二 乘 解 ，SVD 求解 解决 了 舍 入 
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(3) 
(4) 


(5) 


(6) 


| 是 GGT 与 GTG 的 特征 向 量 。diag (wi) 和 


diag(A;)#:4 x 4 的 对 角 阵 ，wi; 和 XA (i = 1,2,3,4) 分 别 是 G 和 G+ 的 奇异 值 ，w1 <w < ws < 
W4» Ay > Az > 43 >A, 


2.1 基于 SVD 方 法 的 均一 方差 定位 误差 模型 


用 a 表示 定位 误差 ， 定 位 误差 协 方差 可 表示 为 : 
Cov(é,) = GtE(Epe," (Gt) = Vdiag(Ai)UTE(epep" )Udiag(wi)VT 


(7) 
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由 式 7) 可 得 定位 误差 的 协 方差 矩阵 ; 


E(e,) = 0 


Cov(é,) = 0a027TNxN 
E(e_.%e,) = 0,1 +j 
Py 69? 为 卫星 的 测 距 误 方差 。 


Cov(er) = Vdiag(A;?)V" oo? 


D 星 的 测 距 误差 服从 独立 、 零 均值 并 具有 相同 方差 的 正 态 分 布 时, 测 距 误差 模型 特性 如 


(8) 


(9) 


该 矩阵 为 一 实 对 称 和 矩阵 ， 可 表示 为 Hao2，H = (GTG)-: 为 几何 精度 因子 矩阵 。 因 此 ， 
H = Vdiag(12)VT (10) 
V 的 四 个 列 向 量 (Vi Vz Vs V4) 是 矩阵 H 的 特征 向 量 ，1 2 > AQ? > Ag? > 44 是 H 的 四 个 特 
征 值 。 定位 误差 的 分 布 是 一 个 四 维 椭 球 ，(Vi V VV4) 定 义 了 误差 椭 球 四 个 主轴 的 方向 , Vi 是 
主轴 ，ooX1 o042 oo13 oo4 是 椭 球 的 半 轴 长 。 
2.2 基于 SVD 方 法 的 非 均 一 方差 定位 误差 模型 
当 卫 星 的 测 距 误 差 无 法 满足 具有 相同 方差 的 正 态 分 布 的 条 件 时 , 用 om 表示 第 n 颗 卫星 的 
测 距 误 差 方 差 ， 测 距 误 差 模 型 特性 如 下 : 
E(ep)= 0 
Cov(é,) = diag(o,7) 7 = 1,2, ..., N (11) 
E(P) = 0,i #j 
由 式 (7) 可 得 定位 误差 的 协 方差 矩阵 : 
Cov(é,) = Vdiag(A,)U'diag(o,7)Udiag(w,)V', i= 1,2,3,4; n = 1,2,...,N (12) 
可 以 看 出 ， 在 卫星 测 距 误差 非 均 一 的 情况 下 ， 定 位 误差 模型 更 加 复杂 
假设 n 颗 卫星 中 测 距 误差 方差 最 小 的 为 00?，o%? > ao2， 将 式 〈12) 中 的 定位 误差 协 方 
差 矩 阵 展开 可 得 总 的 定位 误差 协 方差 矩阵 Cov(s,)7ro0; 


Cov(sr)7 = 
Uni] [Uni]" 
unz | |u , 
Vdiag(4;:?)VT oo + XA Vdiag(A;) ns i diag(A;)V" (on? — 09") 
Una] LUna 


(13) 

Cov(s)7r 由 同方 差 性 定位 误差 与 测 距 误差 方差 大 于 om? 的 测 距 误差 方差 增 量 部 分 
On nace sees Gaia 量 总 和 两 部 分 组 成 。Cov(sr)7 不 仅 与 V, Ay? 
ao2 有关， 还 与 测 距 误差 方差 增 量 om%2? 一 oo2 和 U 和 矩阵 中 第 mn 行 向 量 有 关 。 

从 式 (13) 可 以 看 出 ，Cov(sr)7 是 一 个 实 对 称 窍 阵 。 用 Vr 的 四 个 列 向 量 表 示 Cov(s)7 的 
RIEA rm E72”E73 éra 表示 Cov(e1)zr 的 特征 值 。 根 据 和 矩阵 特征 值 交织 定理 ，&zx1 与 
Ai? 00 ?关系 如 下 : 

ri > W o? > ër” adi "0" > trs 21°00 > Ep” > Ag? 00? 
fiér = Liter Ay? do? + DN Aen” (Gn? — 00°) (14) 


2 
(Dh li Uni Vri)? + (LE Ai Uni Vy) + OL Ai UniVei)? + OE Ai UniVoi)? 


Znx = Cy Ai UniVxi ) /4zn Zny = (Xi di UniVyi ) [Aan (15) 


Znz = (Xi=1iUniVzi) /zn Zip = CPhidiuniVpi) /Azn 
定位 误差 协 方差 矩阵 可 用 特征 值 与 特征 向 量 表示 为 : 
Cov(er)r = Vr diag(Eqi” Vr" (16) 
由 此 可 知 卫 星 测 距 误 差 非 均一 时 ， 定 位 误差 的 分 布 是 一 个 四 维 椭 球 ，Vri Vrz Vra Vira ie 
义 椭 球 的 四 个 主轴 ，Vri 是 主轴 ，éri Ere Erg E74 是 椭 球 的 半 轴 长 。 该 椭 球 相 比 测 距 误 差 均 一 
假设 下 的 定位 误差 椭 球 ， 主轴 袁 一 特定 方 癌 趋 近 , 该 方向 与 星座 分 布 、 测 距 误 差 方 差 增 大 卫 
星 的 位 置 均 有 关 "”。 误 差 椭 球 轴 长 因 测 距 误差 方差 增 量 而 增加 ， 主 轴 增 加 的 最 多 。 当 o,,? 远 
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大 于 oo? 时， Er 将 远大 于 44”go?， 误 差 椭 球 主要 


沿 st 轴 


方差 增 量 将 大 大 拉 长 定位 误差 椭 球 "。 


3 基于 SVD 方法 的 定位 结果 乡 


本 文 利用 北京 站 与 香港 站 的 实测 数据 对 基于 SVD 分 解 的 定位 误差 模型 进行 验证 。 
响 时 的 数据 进行 定位 解 算 , 此 时 卫星 测 昌 


星 信 号 受到 多 径 干 扰 和 电离 层 闪 烁 影 ! 


分 析 


根据 定位 误差 模型 定位 误差 椭 球 3 


= 轴 将 明显 偏转 拉 长 。 


本 次 研究 首先 对 卫星 伪 距 进行 拟 合 去 趋势 处 理 得 到 卫星 测 
对 受到 多 径 干 扰 与 电离 层 内 烁 事件 影响 的 卫星 进行 定位 解 算 ， 
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方向 被 拉 长 。 所 有 卫星 的 测 距 误差 


选取 卫 


距 误差 显著 增 大 ， 


(D ， 得 到 测 距 误 差 均 一 假设 下 定位 误差 椭 球 主轴 的 特征 向 量 与 特征 值 。 


Al (16) ， 得 到 测 距 误 差 非 均一 情 


3. 1 定位 解 算 参 数 
(1) 地 点 坐标 : 

选择 IGS 站 中 北京 站 BJFS 与 香 
位 误差 的 影响 。 


况 下 实际 定位 误差 椭 球 主旨 
比 定位 误差 椭 球 的 特征 向 量 特征 值 与 实际 定位 结 


HR, Suk 


港 站 HKSL 两 个 地 点 ， 


113°55 40'E, 
(2) 卫星 位 置 与 伪 距 测量 值 


利用 2014 年 4 月 25、26、27 日 ， 
年 9 月 22、23、24 日 ，2014 年 12 月 21、22、23 日 


实测 数据 进行 定位 解 算 。 伪 距 测量 


香港 站 处 于 低 纬 地 区 


2014 年 3 月 20、21、22 


误差 。 接 下 来 利用 SVD 方 法 
得 到 实际 定位 


的 特征 向 量 与 特征 值 
定位 误差 椭 球 模型 的 正确 性 。 


二 果 。 利 用 公式 
利用 公式 (12) 
。 通 过 对 


分 别 考察 多 


REM. 


EA Div ALI E A Jae TA PRT RE 
北京 站 的 坐标 为 39"5450N，116"23 30E， 香 港 站 的 坐标 为 2222 19 N， 
， 电 离 层 闪烁 事件 


日 ，2014 年 6 月 20、21、22 日 ，2014 
，2014 年 2 月 27 日 ，2014 年 3 月 10 日 的 GPS 卫星 


值 与 卫星 位 置 由 IGS 网 站 提供 的 采样 率 为 30s 观 测 文 件 与 导 


航 电 文 获得 。 
3. 2 测 距 误差 提取 方法 


伪 距 测量 值 包含 各 种 误差 .按照 来 源 不 同 大 致 可 分 为 三 个 方面 , 与 卫星 有 关 的 误差 , 与 


信号 传播 有 关 的 误差 及 与 接收 机 有 关 的 误差 ,测量 误差 的 两 个 重要 
根据 变化 快慢 这 个 特征 ， 上 述 各 种 误差 可 以 分 为 仿 差 和 噪声 两 类 
变化 较 慢 的 误差 。 偏差 量 的 大 小 具有 一 定 稳 定性 ， 可 
类 误差 ， 具 有 统计 特征 


差 相反 ， 品 声 是 变化 很 快 的 一 
观测 文件 中 记录 了 卫星 的 伪 E 


o 偏差 指 


直接 测量 或 者 月 


2 


特征 是 其 大 小 和 变化 快慢 ， 


的 是 在 一 段 时 间 内 


数学 模型 来 预测 。 
县 体 值 很 难 被 预 讽 


与 偏 
l. 


E 观 测 值 ， 


的 残 差 项 即 是 测 距 误差 en。 


星 到 接收 机 的 几何 距离 和 伪 距 测量 偏差 。 对 


本 次 研究 利用 多 项 式 拟 合 对 卫星 伪 距 随时 间 变 化 


卫星 伪 距 随时 间 变 化 由 


可 知 卫星 伪 距 变化 存在 一 定 趋势 , 趋势 中 包含 卫 


3.2.1 多 径 效 应 对 测 距 误差 的 影响 


多 径 效应 指 的 是 接收 机 天 线 接收 到 GPS 信 


又 可 能 是 经 过 一 


严重 时 还 可 以 导致 码 相 位 失 锁 和 卫星 信号 失踪 。 


现象 。 即 卫星 仰角 越 小 ， 
大 。 


图 1 为 2014 年 4 月 26 日 北京 站 测 得 的 PRN14、15、 


误差 标准 差 o% 随 仰角 变化 图 。 设置 卫星 仰角 截止 角 为 10 。 从 图 1 可 以 看 出 ， 随 仰角 变 小 ， 卫 


号 经 周围 地 物 反射 后 的 信号 ， 而 每 个 
次 或 多 次 反射 后 到 达 天 线 的 。 多 径 效应 会 降低 载波 相位 以 及 伪 距 的 测量 精度 ， 


0. 2 


星 信号 受 多 径 效应 影响 ，anl 
受到 多 径 效应 影响 。 


5m 变 大 到 2m。 


般 来 说 ， 
多 径 效 应 发 生 概率 越 大 ， 伪 距 测量 ; 


18、21、 


低 仰 角 卫 星 信 
EERIK, ME 


可 通过 仰角 大 小 判断 卫星 信 


22、24 六 颗 卫 星 24 小 时 内 测 距 


线 进行 去 趋势 处 理 , 得 到 


| 线 进行 拟 合 去 趋势 处 理 , 用 均 方 根 和 
确定 系数 R 确 定 多 项 式 系 数 ， 实 现 去 趋势 后 测 距 误差 最 准确 。 


反射 信号 


号 更 容易 发 生 多 径 
E 误 差 越 大 ，o 越 


a 


号 测 距 误 差 是 否 
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N 
a 


测 距 误差 标准 差 /m 


BJFS 2014/04/26 


10 20 40 0 60 a 9 
卫星 仰角 /* 

图 1. 测 距 误 差 标准 差 随 仰角 变化 图 

Figl. on elevation-varying graph 


3.2.2 电离 层 内 烁 对 测 距 误差 的 影响 


电离 层 内 烁 导致 电离 层 " 
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P 传 播 的 无 线 电信 号 快速 起 伏 , 严重 时 会 引起 信号 的 中 断 , 是 影 


响 空间 无 线 电 系统 性 能 的 重要 因素 。 低 纬 的 亦 道 异 常 区 古 全 球 范 围 内 电离 层 内 烁 出 现 最 频繁、 


影响 最 严重 的 地 区 之 一 ”。 


ROTI 是 电离 层 TEC 变化 率 指数 ，ROTI RAKE 


电离 层 闪 烁 可 以 引起 用 户 接收 机 测量 误差 增 大 ”。 
闪烁 现象 越 严重 。 图 2 给 出 了 2014 


年 3 月 香港 站 受 电离 层 闪 烁 影响 卫星 的 测 距 误差 标准 差 mm 随 ROTI 的 变化 情况 。 由 图 2 可 见 
on 在 0. 2m 至 3m 间 变 化 ， 且 随 着 ROTI 值 的 增 大 ，om 呈 增 大 趋势 。 可 通过 ROTI 值 判断 卫星 


言 号 测 距 误差 是 否 受 到 电离 层 闪烁 影 响 。 


3.3 实际 定位 结果 分 析 


测 距 误差 标准 差 /m 


HKSL 2014/03 


2 3 4 5 6 
电离 层 TEC 变 化 率 指数 /TECU 


图 2。 测 距 误差 标准 差 随 ROTI 变化 图 
Fig.2 on ROTI-varying graph 


3.3.1 多 径 效应 影响 下 的 定位 误差 分 析 


在 实测 数据 中 , 卫星 星座 分 布 随时 间 变 化 , 选择 在 一 定时 间 段 内 星座 分 布 变化 较 小 的 数 


据 对 定位 误差 模型 进行 验证 。2014 年 4 月 25 日 、26 日 、27 日 三 日 中 ， 每 日 0 点 20 分 到 0 点 35 分 


〈 历 元 40 一 70) ， 北 京 站 可 视 卫 星 的 方位 角 与 仰角 变化 范围 小 于 5 ， 星 座 分 布 基 本 一 致 。 设 


置 卫星 仰角 截止 角 为 10 。 卫 星 仰角 在 10 到 30 之 间 时 ， 卫 星 受 多 径 效应 


影响 。0 点 20 分 一 0 


点 35 分 内 视野 中 共有 PRN14、15、18、21、22、24 六 颗 卫 星 ， 有 且 仅 有 PRN14 一 颗 卫 星 受 到 多 
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径 效 应 影响 。 图 3 给 出 了 2014 年 4 月 26 日 视野 内 每 颗 卫星 全 天 的 or 随时 间 变 化 曲线 ， 可 见 实 
际 定 位 中 ， 卫 星 测 距 误差 是 非 均一 的 。 当 卫星 受到 多 径 干 扰 时 ，o, 显 车 增 大 。 


2 


BJFS 2014/04/26 


PRN14 


1 
2 
1 
2 
1 
2 
1 
2 


测 距 误差 标准 差 /m 


1 
E 
1 
0 


2014 年 4 月 26 日 ] 


12 
时 间 /h 


16 2 24 


[2 FU Ze a HE ee BENT ZEE ihk 


Fig3. o, time-varying curve on April 26, 2014 


在 站 心 坐 标 系 ENU) 下 进行 定位 解 算 ， 坐 标 系 中 E-N 方 向 为 水 平方 向 。 图 4 为 2014 年 4 


月 25 日 、26 日 、27 日 北京 站 定位 结果 在 水 平方 向 上 投影 的 分 布 图 。 


2014/04/25 0:20~0:35(UT) 


2014/04/26 0:20~0:35(UT) 


BJFS 


North/m 


图 4. 


20144 


测 距 误差 均一 假设 下 定位 


误差 椭 球 主轴 方向 
一 一 ”实际 定位 误差 椭 球 主轴 方向 
@ ”每 日 24h 定 位 结果 
e 每 日 0:20"0:35 定 位 结果 
k ”未 受到 多 径 干 扰 的 卫星 
k ”受到 多 径 干 扰 的 卫星 
FE4 月 25 日 、26 日 、27 日 北京 站 定位 结果 分 布 图 


Fig.4 Positioning of BJFS on April 25, 2014, April 26, 2014, April 27, 2014 


将 0 点 20 分 一 0 点 35 分 内 


的 30 个 水 平方 向 定位 结果 用 红 点 在 图 中 标 出 。 全 天 24 小 时 内 2880 


个 定位 结果 用 蓝 点 在 图 中 标 出 。 图 中 卫星 仰角 60、30、0 分 别 对 应 半径 为 2n、4m、6m 的 圆 ， 
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卫星 方位 角 与 定位 误差 方位 对 应 。 将 受到 多 径 干 扰 的 PRN14 号 卫星 15 分 钟 内 的 方位 角 与 仰角 
均值 用 红色 星 形 标 出 ， 其 他 星 均 由 黄色 星 形 表示 。 

30 个 定位 结果 对 应 30 个 几何 和 矩阵 G6。 对 6 甜 阵 进行 SVD 分 解 ， 得 到 每 一 历 元 测 距 误差 均一 
假设 下 定位 误差 椭 球 主轴 的 特征 向 量 WW 和 特征 值 4 。 将 芒 均 值 廊 在 水 平方 向 上 的 投影 用 黄色 
线段 画 出 。 利 用 实测 数据 提取 的 6) 结合 上 文中 式 (12)〉 MR (16) ， 可 求 得 实际 定位 误差 
椭 球 主轴 的 特征 向 量 Vzj 与 特征 值 2x1， 用 红色 线段 表示 向 量 均值 Vi 在 水 平方 向 上 的 投影 ， 
并 用 三 角形 在 红线 上 标记 出 三 倍 特征 值 均值 E71 在 水 平方 向 上 的 投影 值 Er o 

图 5 给 出 了 2014 年 春分 、 夏 至 、 秋 分 、 冬 至 及 其 前 后 两 天 北京 站 水 平方 向 的 定位 分 布 。 
图 5(a)5(b)5(c)5(d) 4) Hil H20144F3 H 20H ~22 4 04: 10~04:25, 2014476 20H ~22 
A4% A02:30~02:45, 2014469 A 22H ~24H 4H01:35~01:50. 20144712 A 21H ~23 H fi 
日 04:45 一 05:00 的 定位 结果 构成 。 图 中 卫星 仰角 60 、30 、0 分 别 对 应 半径 为 8/3m、16/3m、 
gm 的 圆 ， 卫 星 方 位 角 与 定位 误差 方位 对 应 。 
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(a) 2014/03/20~2014/03/22 04:10~04:25(UT) (b) 2014/06/20~2014/06/22 02:30~02:45(UT) 
8 8 
6 6 
4 4 
2 2 
E £ 
£ 0 £ 0 
o o 
= 三 
-2 -2 
-4 -4 
-6 -6 
-8 -8 
-8 -6 -4 -2 0 2 4 6 8 -8 -6 -4 -2 0 2 4 6 8 
East/m East/m 
(c) 2014/09/22 ~2014/09/24 01:35~01:50(UT) (d) 2014/12/21~2014/12/23 04:45~05:00(UT) 
8 8 
6 6 
4 4 
2 2 
£ £ 
g ° € 0 
9 o 
2 2 
-2 -2 
-4 -4 
-6 -6 
-8 -8 
-8 -6 -4 -2 0 2 4 6 8 -8 -6 -4 -2 0 2 4 6 8 
East/m East/m 
图 5.。 ”春分 、 夏 至 、 秋 分、 冬至 北京 站 定位 结果 
Fig5. Positioning of BJFS on 2014/3/20°22, 2014/6/20°22, 2014/9/22°24, 2014/12/21°23 


表 1 中 给 出 了 2014 年 3 月 21 日 04:10 一 04:25、2014 年 4 月 26 日 0: 20~0: 35, 20144F6 


21 


日 02:30~02:45、2014 年 9 月 23 日 01:35~01:50、2014 年 12 月 22 日 04:45~05:00 每 日 30 个 历 
元 内 的 儿 何 精度 因子 G6D0P 均 值 , 测 距 误差 方差 增 量 均值 (on? 一 


小 值 )， 测 距 误差 均一 假设 下 定位 误差 椭 球 主轴 


BY /二 J 
AEST © 


ENE TR Ze MER EAE RE EL ACP 


表 1 北京 站 定位 结果 分 析 
Tab.1 Analysis of Positioning of BJFS 


00") 《oo 取 视 野 内 卫星 om 最 


轴 长 oo 均值 的 水 平 投影 值 co 和 和 实际 
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几何 精度 因子 。 测 距 误 差 方 差 ” 测 距 误 差 均 一 假设 下 定位 误差 。 实际 定位 误差 椭 球 主轴 半 轴 


作 期 刊 


均值 /m 增 量 均值 /m ” 椭 球 主轴 半 轴 长 均值 投影 值 /m 长 均值 投影 值 /m 
03/21 04:10~04:25 3. 594 2. 587 0. 470 2. 078 
04/26 0:20 ~ 0:35 3. 598 2. 655 0. 495 2. 092 
06/21 02:30~02:45 2. 707 4.811 0. 396 1. 937 
09/23 01:35~01:50 2. 654 3. 495 0.417 1, 822 
12/22 04:45~05:00 2. 853 2. 875 0. 447 1. 986 


结合 图 、 表 可 以 看 出 ,实际 定位 结果 与 非 均一 方差 定位 误差 模型 基本 吻合 ,在 轴 方 向 上 ， 
实际 定位 误差 椭 球 主轴 方向 相 比 测 距 误 差 均 一 假设 下 的 误差 椭 球 主轴 方向 明显 旋转 , 旋转 角 
度 与 星座 分 布 、o% 增 大 卫星 位 置 有 关 ， 所 选 时 间 段 内 的 定位 点 基本 分 布 在 实际 定位 误差 椭 
球 主轴 方向 周围 。 对 于 轴 长 , 实际 定位 误差 椭 球 轴 长 相 比 测 距 误 差 均一 假设 下 的 误差 椭 球 轴 
长 明显 增 大 ， 增 大 量 与 星座 分 布 和 测 距 误 差 方 差 增 量 有 关 ， 实 际 定位 结果 分 布 在 3571 范围 
内 ， 符 合 中 心 极限 定理 。 
3.3.2 电离 层 闪 烁 影响 下 的 定位 误差 分 析 
通过 ROTI 值 可 以 判断 卫星 信号 是 否 受 到 电离 层 内 烁 影响 。 根 据 计 算 结果 得 知 ，2014 年 2 
月 27 日 16:00 一 17:00PRN7 号 、9 号 卫星 与 2014 年 3 月 10 日 18:00 一 19:00PRN7 号 卫星 均 受 到 较 强 
的 电离 层 内 烁 事件 影响 。 

其 中 ,2014 年 2 月 27 日 16 时 25 分 一 16 时 40 分 卫星 星座 变化 较 小 , 此 时 视野 内 共有 PRN1、7、 
9、11、19、28、32 七 颗 卫 星 ， 其 中 PRN19、32 号 卫星 受到 多 径 干 扰 ，PRN7 号 、9 号 卫星 受到 
BAJE NERI o 

KeA HY TREMPE Yo, BAN Ta) ee HZ. FS 4 MS 82 Be Ey 
层 闪 烁 影响 时 ，o% 明 显 增 大 。 
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图 6，2014 年 2 月 27 日 香港 站 卫星 测 距 误 差 标准 差 随 时 间 变 化 曲线 


Fig6. o, time-varying curve on Feb 27, 20 
图 7 为 2014 年 2 月 27 日 香港 站 水 平方 向 的 定位 结果 分 布 图 。 蓝 点 标记 全 天 定位 结果 。 红 点 

标记 16 时 25 分 一 16 时 40 分 内 的 定位 结果 。 用 红色 星 形 标记 受到 多 径 干 扰 的 PRN19 号 卫星 与 
PRN32 号 卫星 。 用 绿色 星 形 标记 受到 电离 层 内 烁 影响 的 PRN7 号 卫星 与 PRN9 号 卫星 。 其 他 星 均 
昌 黄 色 星 形 标记 。 

经 计算 得 ， 嫩 7 =2. 3149m, oA, =0. 4119m。 黄 色 线 段 与 红色 线段 分 别 代表 测 距 误差 均 
一 假设 下 误差 椭 球 主轴 方向 与 实际 定位 误差 椭 球 主轴 方向 在 水 平方 向 上 的 投影 ,用 三 角形 在 
实际 定位 误差 椭 球 主轴 方向 上 标记 出 三 倍 厨 7 。 
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图 8 为 2014 年 3 月 10 日 香港 站 水 平方 向 的 定位 结果 分 布 图 .17 时 15 分 一 17 时 30 分 卫星 星座 
变化 较 小 ， 用 红 点 标记 这 段 时 间 内 的 定位 结果 。 这 段 时 间 视 野 内 有 PRN1、4、7、8、9、11、 
17、20、28、32 十 颗 卫 星 。PRN7、11、32 三 颗 卫 星 受到 多 径 干 扰 ， 其 中 PRN7 号 卫星 在 这 段 时 
间 内 既 受 到 多 径 干扰 又 受到 电离 层 内 烁 影响 ， 用 黑色 星 形 标记 。 经 计算 得 ，571 =2. 1574m, 
gA, =0.4072m。 


2014/02/27 16:25~16:40(UT) 2014/03/10 17:15~17:30(UT) 


8 HKSL 


North/m 
North/m 


East/m 


图 7.、 ”2014 年 2 月 27 日 香港 站 定位 结果 图 8.， ”2014 年 3 月 10 日 香港 站 定位 结果 
Fig.7 Positioning of HKSL on Feb 27, 2014 Fig.8 Positioning of HKSL on Mar 10, 2014 
结合 定位 分 布 图 与 计算 结果 可 知 , 实际 定位 结果 在 轴 疝 、 轴 长 两 方面 均 与 非 均 一 方差 定 
位 误差 模型 相 吻 合 。 且 多 颗 卫 星 测 距 误差 标准 差 均 显著 增 大 时 ,误差 椭 球 主轴 拉 长 更 为 明显 。 


4 结论 


本 文 对 GNSS 卫 星 测 距 误 差 非 均 一 时 的 定位 误差 特性 进行 研究 ， 利 用 SVD 分 解 方 法 给 出 了 
非 均 一 方差 定位 误差 模型 。 模 型 指出 ， 当 测 距 误差 非 均一 时 ， 定 位 误差 在 空间 中 呈 四 维 椭 球 
分 布 ， 椭 球 的 轴 向 和 轴 长 可 由 对 定位 误差 协 方差 矩阵 进行 SVD 分 解 获 得 。 与 测 距 误差 均一 假 
设 下 的 定位 误差 椭 球 相 比 ， 测 距 误差 非 均一 时 定位 误差 椭 球 因 测 距 误 差 方差 增 量 偏转 拉 长 。 
当 有 卫星 测 距 误差 方差 显著 增 大 时 , 主轴 偏转 明显 ， 椭 球 主要 沿 主 轴 被 拉 长 。 本 文 利用 北京 
站 与 香港 站 的 实测 数据 对 模型 进行 检验 。 首先 对 卫星 伪 距 进行 拟 合 去 趋势 处 理 , 提取 出 卫星 
测 距 误差 ; 接着 对 受 多 径 干 扰 和 电离 层 闪 烁 影响 的 卫星 进行 定位 解 算 分 析 。 结 果 表 明 : 卫星 
审 号 受到 多 径 干 扰 或 电离 层 闪 烁 影响 时 , 测 距 误差 标准 差 显 著 增 大 。 此 时 实际 定位 误差 椭 球 
相 较 测 距 误差 均一 假设 下 的 定位 误差 椭 球 ,主轴 方向 明显 旋转 ,主轴 明显 拉 长 。 实 际 定位 结 
果 与 非 均一 定位 差 椭 球 主轴 方向 一 致 ， 分 布 范 围 与 椭 球 主轴 半 轴 长 吻合 。 可 见 ， 非 均一 方差 
定位 误差 模型 可 以 准确 描述 几何 精度 因子 、 测 距 误差 与 定位 误差 之 间 的 关系 , 为 进一步 研究 
定位 精度 提供 依据 。 


a 
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A test of GNSS positioning error model of 


non-uniform pseudorange error 
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Abstract: In the global navigation satellite system (GNSS), an accurate positioning error 
model is important for understanding how the system works . When the pseudorange error of each 
satellite has an uncorrelated normal distribution with zero-mean and common variance , the 
positioning error is equal to the standard deviation of pseudorange error multiplied by the 
geometric precision. However, in practical applications, the pseudorange errors may differ in 
variance due to the different signal propagation paths. Especially when the signals are affected by 
multipath effects or ionospheric scintillation, the pseudorange errors will be significantly different 
from one another. This paper introduces a theoretical positioning error model of non-uniform 
pseudorange error based on singular value decomposition (SVD). Using GNSS data at Beijing 
station and Hong Kong station to analyze the positioning errors affected by multipath effects and 
ionospheric scintillation. The statistical characteristics of positioning error are consistent with the 
results given by the positioning error model of non-uniform pseudorange error, proving the 
correctness and validity of the model. 

Keywords: Positioning error; Geometric precision; Pseudorange error; Multipath effect; 


Ionospheric scintillation 


